Superheavy magic structures and relativistic pseudo-spin symmetry 
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We have explored the occurrence of the spherical shell closures for superheavy nuclei in the framework of the 
relativistic Hartree-Fock-Bogoliubov (RHFB) theory and the shell effects in terms of two-nucleon gaps S^ij,). 
Although the results depend slightly on the forces used, the general set of magic numbers beyond 208 Pb are 
predicted as Z = 120, 138 for protons and N = 172, 184, 228 and 258 for neutrons, respectively. Specifically the 
RHFB calculations favor the nuclide 304 120 as the next spherical doubly magic one beyond 208 Pb. As indicated 
by the self-consistent calculations, the occurrences of superheavy shell closures are essentially related with the 
violation and restoration of relativistic pseudo-spin symmetry. Further analysis imply that the shell effects are 
sensitive to both the values of scalar and nonrelativistic effective masses, which togetherwith the violation of 
the relativistic symmetry may interpret the origin of model deviations. 
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For a fairly long period, it remains a challenging issue in 
nuclear physics to explore the mass and charge limit of realis- 
tic nuclei, i.e., the superheavy elements (SHEs) with Z > 104 
and the so-called stability island of superheavy nuclei (SHN). 
Far from being simply large clusters of nucleons, these fas- 
cinating species owe their very existence to subtle contribu- 
tions to the nuclear binding energy [l|j. The theoretical stud- 
ies, including the macroscopic-microscopic models and 
self-consistent mean field approaches HE), have provided a 
large amount of valuable information for the exploration of 
SHN. Many efforts have been devoted to finding the next sta- 
bility island, i.e., the magic SHN with the enhanced stabil- 
ity promised by the shell effects. Although the consensus 
is somehow reached on the spherical neutron shell closure 
N - 184, other possible neutron magic shells such asN - 172 
and the next spherical proton shell probably lying in the range 
of 1 14 < Z < 126 are still being under discussions J6H7J. 

Experimentally, the discoveries of new elements up to Z = 
118 have been reported in Refs. [8] [9]. However, these SHN 
are not located on the stability island as originally sought. 
In addition the enlarged survival probabilities with increasing 
proton number of SHEs from Z = 1 14 to 1 18 seem to indicate 
the enhanced shell effects with increasing Z and therefore a 
possible proton magic shell may emerge beyond Z = 1 18 IflOl . 

On the other hand, the extrapolation towards the super- 
heavy region also challenges the predictive reliability of nu- 
clear structure models. The macroscopic-microscopic meth- 
ods, although generally successful in describing the nuclear 
binding, require preconceived knowledge of the expected den- 
sities and single-particle (s.p.) potentials J2), which fades 
away when stepping into new regions where stronger polar- 
ization effects and more complicated functional forms of the 
densities may occur J4|. It is well known that the SHN are 
mainly stabilized by the shell effects, and the self-consistency 
of the theory is therefore of special significance for a reliable 
description of SHN. In fact, the shell effects are essentially 
related with the nuclear spin-orbit(SO) splittings and relevant 
pseudo-spin symmetry restoration ifTTI as well as the shape of 
the nucleus. Whereas the variations of SO splitting are mainly 



determined by both SO potential concerned with the derivative 
of the average potential and the tensor force originated from 
the exchange of pion lfl2l . This leads to a challenging and 
ambitious task to provide a unified self-consistent description 
of both spin-orbit coupling and tensor effects, which can be 
achieved on the level of the relativistic Hartree-Fock (RHF) 
approach under the meson exchange picture of nuclear forces. 

In the past years, the covariant density functional (CDF) 
theories lfL3l [14 1 have achieved great success in the self- 
consistent description of nuclear systems, including the sta- 
ble and weakly bound nuclei and neutron stars. One of exist- 
ing approaches is the density-dependent relativistic Hartree- 
Fock (DDRHF) theory lfT31 [TBI and its extension — the rela- 
tivistic Hartree-Fock-Bogoliubov (RHFB) theory lfl7l . Within 
DDRHF, the Lorentz covariant structure is fully kept, thus en- 
suring the natural treatment of nuclear spin-orbit couplings. 
The inclusion of Fock terms in DDRHF not only brings sub- 
stantial improvements in self-consistent descriptions of nu- 
clear structure properties [18-20] and neutron star physics 
ET1 [22], but also provides the natural origin of the nuclear 
tensor force, e.g., the pion exchange and p-tensor couplings 

HEED- 

As an important residual interaction beyond the p-h chan- 
nel, the pairing correlations not only contribute to the stability 
of SHN but also play a major role in determining the nuclear 
deformation Il23ll . which requires a unified self-consistent de- 
scription of the mean field, pairing field and nuclear defor- 
mation. In terms of Bogoliubov quasi-particles, such a goal 
is partially achieved by the RHFB theory which unifies the 
description of RHF mean field and pairing correlations. Al- 
though a complete RHFB study would require to explore also 
the shape degrees of freedom [24] — a very time-consuming 
task — one can already discuss the spherical magicity of SHN 
in the spherical RHFB calculations. 

In this work we investigate the superheavy nuclides cover- 
ing Z = 110 ~ 140. In the pairing channel, the finite-range 
Gogny force D1S l25l renormalized by a strength factor / 
is adopted as effective pairing interaction. The introduction of 
the factor / is due to the fact that the pairing related quantities, 
such as the odd-even mass differences and the moments of in- 
ertia, are systematically overestimated in the RHFB calcula- 
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tions of heavy nuclei with the original Gogny pairing force 
11261 . To eliminate such systematical deviations, the strength 
factor is reduced as / = 0.9 by reproducing the odd-even mass 
differences of odd Pb isotopes. Concerning the RHF mean 
field, the adopted effective Lagrangians are PKA1 lfl6l and 
PKO series lfT5l[T8l . For comparison the relativistic Hartree- 
Bogoliubov (RHB) calculations HU with PKDD (27) and 
DD-ME2 l28ll are also presented in this work. The integro- 
differential RHFB equations are solved on a Dirac Woods- 
Saxon basis [29 1 with a radial cutoff R — 28 fm. The numbers 
of positive and negative energy states in the basis expansion 
for each s.p. angular momentum (/, j) are chosen to be 44 and 
12, respectively. 

To identify the magic shells, different observables can be 
chosen, such as nucleon separation energies, pairing gap en- 
ergy, s.p. spectra, etc. The magic shells in SHN might not 
be as well-marked as in the ordinary nuclei. Here, we con- 
centrate on two types of observables. One is the so-called 
two-nucleon gaps [6 2p (proton) and 6 2 „ (neutron)], i.e., the 
difference of two-nucleon separation energies of neighboring 
isotopes or isotones, which provides an efficient evaluation of 
the shell effects, 

S 2p (N,Z) = S 2p (N,Z)-S 2p (N,Z + 2), (la) 
5 2n (N,Z) = S 2n (N,Z)-S 2n (N + 2,Z). (lb) 

Obviously, the peak values of the two-nucleon gaps are es- 
sentially determined by the sudden jump of the two-nucleon 
separation energies, which can be taken as the evidence of the 
magic shell occurrence. 

Fig. [T]presents the two-proton (left panels) and two-neutron 
(right panels) gaps for the even-even isotopes of Z = 110 ~ 
140, calculated with the selected effective Lagrangians. The 
thick stepped lines show the two-proton drip lines identified 
by two-proton separation energies. Nuclei that are stable with 
respect to /3-decay and fission in our calculations are em- 
phasized by filled stars and circles, respectively. The /? (or 
fission)-stability line is defined by the maximum binding en- 
ergy per nucleon for a given A (or Z), which corresponds to 
the minimum Q-value of /3-decay (or fission) [31 1. The dashed 
lines denote the empirical /^-stability lines described by the 
formula Z = A/(1.98 + 0.0155A 2/3 ) [30). The open thick 
squares lying at the lower-left corner of each plot represent 
the even-even SHN observed or declared to be observed ex- 
perimentally. It is clearly seen that the experimentally discov- 
ered elements coincide remarkably well with the predicted fis- 
sion stability line, especially for PKA1, although the deforma- 
tion may play a significant role in stabilizing the SHN. Strictly 
speaking, only nuclei having negative chemical potentials are 
stable against nucleon emission. However, due to the strongly 
enhanced Coulomb barriers in SHN the two-proton drip lines 
extend by a few more units to the proton-rich side beyond 
negative chemical potentials. In Fig{T] the squares filled with 
thicker color stand for the larger gaps and obviously the peak 
areas indicate the emergence of magic shells. 

From Fig. [T]one can see that PKA1 presents the most dis- 
tinct shell effects, supporting Z = 120, 126, 138 and N - 184, 
258 as the candidates of proton and neutron magic numbers, 
respectively. The other effective Lagrangians also present 



a fairly remarkable proton shell at Z = 120. In addition, 
Z = 132 by PKDD and Z = 138 by both RHFB and RCHB 
calculations are found to be possible proton magic numbers, 
consistent with the predictions in Ref. [5|. Concerning the 
neutron shells, besides N = 184 and 258, PKA1 also presents 
fairly distinct shell structures atN = 172, which is illustrated 
evidently by the other Lagrangians. Fairly distinct shell effects 
atN = 184 and 258 are also found with the other parametriza- 
tions except with PK02. Remarkable shell effects are also 
found atN — 228, although not so strongly as at N = 184 and 
258 by PKA1. Furthermore, the occurrence of neutron shells 
at N = 164 is illustrated by PK02, PKDD and DD-ME2, and 
the latter two demonstrate another occurrence atN = 198. 

Inspecting another observable — the pairing gap energies 
whose zero values accompany the emergence of magic shells, 
one can achieve similar conclusions as above on the emer- 
gence of proton and neutron shells. Combined with the two- 
nucleon gaps, it is found that the proton shell Z = 120 is 
supported evidently by PKA1, and by the other interactions 
to some extent. Hence Z = 120 can be considered as model- 
independent candidate of proton magic numbers. On the other 
hand, the situation for the neutrons is more complex. Al- 
though N - 112 and 228 seem to be generally supported by 
the selected effective Lagrangians, the corresponding shell ef- 
fects are not strong enough to determine the occurrence of 
magic shells. Except for PK02, N = 184 and 258 are also 
generally favored as the candidates of neutron magic numbers. 
Specifically, PKA1 can provide a better description on the nu- 
clear shell structure than the others [ 16 1 and better agreement 
on the fission stability of discovered SHN (see Fig. [TJ, and 
presents pronounced shell effects there. In fact, as indicated 
by Skyrme Hartree-Fock (SHF) investigations E21 N = 184 
is also favored evidently to be a spherical neutron magic num- 
ber and the N = 184 isotones are expected to have spherical 
shapes. 

TABLE I: Bulk properties of symmetric nuclear matter calculated 
with the effective interactions PKA1, PKO series, PKDD and DD- 
ME2: saturation density p (fm -3 ), binding energy per particle Eg/A 
(MeV), incompressibility A* (MeV), asymmetry energy coefficient J 
(MeV), scalar mass M* s and non-relativistic effective mass M* NR 1151 
in units of nucleon mass M. 



Force 


Po 


E B /A 


K 


J 


M* 


M NR 


PKA1 


0.160 


-15.83 


229.96 


36.02 


0.547 


0.681 


PKOl 


0.152 


-16.00 


250.24 


34.37 


0.590 


0.746 


PK02 


0.151 


-16.03 


249.60 


32.49 


0.603 


0.764 


PK03 


0.153 


-16.04 


262.47 


32.98 


0.586 


0.742 


PKDD 


0.150 


-16.27 


262.19 


36.79 


0.571 


0.651 


DD-ME2 


0.152 


-16.11 


250.30 


32.27 


0.572 


0.652 



Nevertheless, from Fig. [T] one can find distinct devia- 
tions among the models in predicting the magic numbers. 
Z = 120 can be considered as a reliable prediction of pro- 
ton magic number and Z = 138 could be another candidate 
with more model dependence. The neutron shells N = 172, 
184, 228 and 258 are common to several models. Other 
shells, e.g., N = 198, appear essentially model dependent. 
Among the present models, one may notice that RHB calcula- 
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FIG. 1 : (Color online) Contour plots (in MeV) of two-proton gaps 62 P (left panels) and two-neutron gaps 62* (right panel) as functions of N and 
Z. The results are extracted from the RHFB calculations with PKA1, PK02 and PK03, and RHB ones with PKDD and DD-ME2. The thick 
stepped lines show the two-proton drip lines. The nuclei stable with respect to /3-decay and fission are marked with filled stars and circles, 
respectively. The long-dashed lines represent the yS-stability line Z = A/( 1.98 + 0.01 55A 2/3 ) [ 30 1 . The open thick squares lying at the lower-left 
corners denote the superheavy nuclei observed or declared to be observed experimentally. 



tions (PKDD and DD-ME2) predict more shell closures than 
the others while PK02 predicts fewest candidates. To inter- 
pret such distinct deviations, Table [I] shows the bulk proper- 
ties of symmetric nuclear matter determined by the present 
sets of Lagrangians. In general the occurrence of superheavy 
magic shells is closely related with both the scalar mass M* s 
and effective mass M* NR 1151 . which essentially determine the 
strength of spin-orbit couplings and level densities, respec- 
tively. Among the present models, the effective Lagrangian 
PK02 predicts the largest values of both masses, which leads 
to relatively weak spin-orbit couplings and high level density 
on the average. As a result there remains little space in the 
spectra for the occurrence of magic shells. On the other hand, 
the RHB models predict more magic shells due to the rela- 
tively small masses. In fact, as seen from Fig. [T] PK02 also 
presents weaker shell effects than the others. For PKA1 the 
situation is different. Although it has a larger effective mass 
M* NR than PKDD or DD-ME2, PKA1 gives a smaller scalar 
mass M* s and shows stronger shell effects than the others. Due 
to these reasons PKA1 does not suffer of the common draw- 
back of the CDF calculations — the so-called artificial shell 
closures ||33l — and it better preserves the relativistic sym- 
metry lfl6l l20l . From this point of view the deviations be- 
tween models may also originate from the relativistic symme- 
try restorations. 

Similarly to the situation in the stable region [ 16], the model 



deviations originating from the relativistic pseudo-spin sym- 
metry conservation can also be found in the s.p. spectra of 
SHN. Taking the doubly magic SHN 304 1 20i8 4 as an example, 
Fig. [2] shows the proton (left panel) and neutron (right panel) 
canonical s.p. spectra provided by selected models. It is found 
that PKA1 provides the most evident magicity at Z = 120 
and N — 184, respectively, although these shell closures are 
much weaker than in ordinary nuclei. For the neutron shell 
N — 184, it is essentially determined by the degeneracy of two 
pseudo-spin partners [2hu/2, 1 713/2} an d {4*i/2> 3^3/2}, respec- 
tively above and below the shell. For the latter, the pseudo- 
spin symmetry is well preserved in all the calculations while 
for the former with high angular momentum the symmetry 
is only properly restored by PKA1 while seriously violated 
by the others, leading to the occurrence of the shell closure 
N = 198. In fact, a similar emergence of shell closure can 
also be found in mid-heavy and heavy regions of the nuclear 
chart. For instance, the proton shell Z = 82 in 208 Pb can be 
also taken as the result of the degeneracies of two pseudo-spin 
partners {2/7/2, 1/29/2} and {3*1/2,2^3/2}. In the CDF calcula- 
tions (except PKA1) there is a clear gap between 2/7/2 and 
I&9/2, i.e., the artificial shell closure Z = 92 11331 , which some- 
what compresses the magic shell Z = 82. 

A similar mechanism can also be found in the formation of 
sub-shell 64 due to the degeneracy of the pseudo-spin partners 
{3ii/2, 2^3/2} and {2c/s/2, Igi/i] [16|. In the CDF calculations 
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FIG. 2: (Color online) Proton (left panel) and neutron (right panel) canonical s.p. spectra of superheavy nuclide 304 120i84. The results are 
extracted from the RHFB calculations with PKO series and PKA1, as compared to the RHB ones with PKDD and DD-ME2. In all cases the 
pairing force is taken as the finite range Gogny force D1S with the strength factor / = 0.9. See the text for details. 



(except with PKA1) the sub-shell 64 is compressed by the vio- 
lation of pseudo-spin symmetry on the partners {2ds/2, igi/i}, 
which induces the so-called artificial shell closure 58 lfl6l[33l . 
From this point of view, the restoration of pseudo-spin sym- 
metry may also play a delicate role for the occurrence of shell 
closures in SHN, and this could be another origin of the devi- 
ations between the models. 

In fact, such tight relation between the relativistic pseudo- 
spin symmetry restoration and the occurrence of magic shells 
is not occasionally existing in the superheavy region. It is in- 
teresting to notice that there exists at least one pair of pseudo- 
spin partners above or below the conventional magic shells 
except N/Z = 2 and 8 and the occurrences of shell closures 
are also favored by the relativistic symmetry restoration to dif- 
ferent extent. Traditionally these magicities are interpreted by 
the strong spin-orbit effects of high- j orbits. While for SHN 
the increased s.p. levels will 'intrude' the space left by the 
spin-orbit splitting of high- j orbits such that the restoration of 
relativistic pseudo-spin symmetry may play more significant 
roles in determining the occurrence of shell closures. The typ- 
ical example is the neutron one at N = 184 demonstrated by 
PKA1 (see the right panel of Fig. [2}. 

Different from N = 184, which is favored by the well- 
preserved pseudo-spin symmetries, the proton shell closure 
Z = 120 emerges with the violation of pseudo-spin symmetry. 
As shown in the left panel of Fig. [2] the proton shell closure 
Z = 120 is determined directly by the large splittings of two 
pseudo-spin partner states, 3p^/2 and 2/5/2, whereas the spin 
doublet 3p above the shell is nearly degenerated, which can be 
qualitatively interpreted by the proton configuration. Below 
the shell Z = 120, the protons filling in the high- /' states will 
be driven from the center of nucleus due to the strong centrifu- 
gal potential. In addition, the enhanced repulsive Coulomb 
field in SHN will also expel the protons to spread over large 



radial distance. Both effects lead to an interior depression of 
the proton distributions and consequently the interior region 
of the mean potential is not flat any more, which may reduce 
(enlarge) the spin (pseudo-spin) orbital splitting, particularly 
for low- j states 3p and 2f. In Ref. [34| it is also pointed out 
that the distinct central depressions on the densities lead to the 
spherical shell gaps at Z = 120 and N - 172 as a direct con- 
sequence of pseudospin symmetry breaking, whereas a flatter 
density profile favors the shell occurrence atN= 184 as well 
as the proton one at Z = 126. In fact not only for SHN, the 
emergence of new shell closure Z/N = 16 and N -32 Il35ll36ll 
can be also related with the violation of pseudo-spin symme- 
try in light exotic nuclei. 

In summary, we have explored the occurrence of spherical 
shell closures for superheavy nuclei and the physics therein 
using the relativistic Hartree-Fock-Bogoliubov (RHFB) the- 
ory with density-dependent meson-nucleon couplings. The 
shell effects are quantified in terms of two-nucleon gaps 52n( P )- 
To our knowledge, this is the first attempt to perform such 
extensive calculations within the RHFB scheme. The results 
indicate that the nuclide 304 120 could be the next spherically 
doubly magic nuclide beyond 208 Pb. It is also found that the 
shell effects in SHN are sensitive to the values of both scalar 
mass and effective mass, which essentially determine the spin- 
orbit effects and level density, respectively. In addition, the 
emergence of shell closures is found to be concerned with the 
violation and restoration of relativistic pseudo-spin symme- 
try. To some extent, one may further conclude that the model 
deviations may originate from both effective masses and rel- 
ativistic symmetry. Experimental measurement of Q a for at 
least one isotope of Z = 120 nucleus would help us to set 
proper constraint in determining the shell effects of SHN and 
further test the reliability of the models as well. One also has 
to admit that for a more extensive exploration one needs to 
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take into account the deformation effects. 

One of the authors (J. J. Li) would like to express his grati- 
tude to Mr. H. F. Zhang for his helpful discussions. This work 
is partly supported by the National Natural Science Founda- 



tion of China under Grant No. 11075066, the Fundamental 
Research Funds for the Central Universities under Contracts 
No. Izujbky-2012-k07, and the Program for New Century Ex- 
cellent Talents in University. 



[1] M. Block et al., Nature 463, 785 (2010). 

[2] P. Moller, J. Nix, W. Myers, and W. Swiatecki, At. Data Nucl. 

Data Tables 59, 185 (1995). 
[3] A. Baran, Z. Lojewski, K. Sieja, and M. Kowal, Phys. Rev. C 

72,044310(2005). 
[4] M. Bender, K. Rutz, P.-G. Reinhard, J. Maruhn, and W. Greiner, 

Phys. Rev. C 60, 034304 (1999). 
[5] W. Zhang, J. Meng, S. Zhang, L. Geng, and H. Toki, Nucl. 

Phys. A 753, 106 (2005). 
[6] R. Gupta, M. Manhas, W. Greiner, et al., J. Phys. G: Nucl. Part. 

Phys 36, 115105 (2009). 
[7] H. Zhang, Y. Gao, N. Wang, J. Li, E. Zhao, G. Royer, et al., 

Phys. Rev. C 85, 014325 (2012). 
[8] Y. T. Oganessian et al., Phys. Rev. C 74, 044602 (2006). 
[9] Y. T. Oganessian et al., Phys. Rev. Lett 104, 142502 (2010). 
[10] G. Adamian, N. Antonenko, V. Sargsyan, and W. Scheid, Nucl. 

Phys. A 834, 345c (2010). 
[11] W. Long, T. Nakatsukasa, H. Sagawa, J. Meng, H. Nakada, and 

Y. Zhang, Phys. Lett. B 680, 428 (2009). 
[12] T. Otsuka, T. Suzuki, R. Fujimoto, H. Grawe, and Y. Akaishi, 

Phys. Rev. Lett. 95, 232502 (2005). 
[13] P.-G. Reinhard, Rep. Prog. Phys 52, 439 (1989). 
[14] J. Meng, H. Toki, S. Zhou, S. Zhang, W. Long, and L. Geng, 

Prog. Part. Nucl. Phys 57, 470 (2006). 
[15] W. H. Long, N. Van Giai, and J. Meng, Phys. Lett. B 640, 150 
(2006). 

[16] W. H. Long, H. Sagawa, N. Van Giai, and J. Meng, Phys. Rev. 

C 76, 034314(2007). 
[17] W. H. Long, P. Ring, N. Van Giai, and J. Meng, Phys. Rev. C 

81, 024308 (2010). 

[18] W. Long, H. Sagawa, J. Meng, and N. Van Giai, Europhys. Lett 

82, 12001 (2008). 

[19] H. Liang, N. Van Giai, and J. Meng, Phys. Rev. Lett. 101, 
122502 (2008). 

[20] W. H. Long, P. Ring, J. Meng, N. Van Giai, and C. Bertulani, 



Phys. Rev. C 81, 031302 (2010). 
[21] B. Sun, W. Long, J. Meng, and U. Lombardo, Phys. Rev. C 78, 
065805 (2008). 

[22] W. H. Long, B. Y. Sun, K. Hagino, and H. Sagawa, Phys. Rev. 

C 85, 025806 (2012). 
[23] S. K. Patra, C.-L. Wu, C. R. Praharaj, and R. K. Gupta, Nucl. 

Phys. A 651, 117(1999). 
[24] J. -P. Ebran, E. Khan, D. Pena Arteaga, and D. Vretenar, Phys. 

Rev. C 83, 064323 (2011). 
[25] J. Berger, M. Girod, and D. Gogny, Nucl. Phys. A 428, 23 

(1984). 

[26] L. J. Wang, B. Y. Sun, J. M. Dong, and W. H. Long, submitted 
to Phys. Rev. C (2013), arXiv: 1302.5817 [nucl-th]. 

[27] W. Long, J. Meng, N. Van Giai, and S. Zhou, Phys. Rev. C 69, 
034319 (2004). 

[28] G. Lalazissis, T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. C 

71,24312(2005). 
[29] S.-G. Zhou, J. Meng, and P. Ring, Phys. Rev. C 68, 034323 

(2003). 

[30] P. Marmier and E. Sheldon, Physics of Nuclei and Particles , 

vol. 1 (Academic Press, New York, 1971). 
[31] C.-L. Wu, M. Guidry, and D. H. Feng, Phys. Lett. B 387, 449 

(1996). 

[32] A. Kruppa, M. Bender, W. Nazarewicz, P. Reinhard, T. Vertse, 
and S. Cwiok, Phys. Rev. C 61, 034313 (2000). 

[33] L.-S. Geng, J. Meng, H. Toki, W.-H. Long, and G. Shen, Chin. 
Phys. Lett 23, 1139 (2006). 

[34] A. Afanasjev and S. Frauendorf, Phys. Rev. C 71, 024308 
(2005). 

[35] A. Ozawa, T. Kobayashi, T. Suzuki, K. Yoshida, and I. Tanihata, 

Phys. Rev. Lett. 84, 5493 (2000). 
[36] R. Kanungo, I. Tanihata, and A. Ozawa, Phys. Lett. B 528, 58 

(2002). 



